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Lecture 6 



Assembly of independent molecules 

So far the description of the interaction between light and matter has been in a very general form, 
where it was assumed merely that the interaction is local, and in the electric dipolar approximation, 
where magnetic dipolar and electric quadrupolar interactions (as well as higher order terms) were 
neglected. The ensemble of molecules has so far no constraints in terms of composition or mutual 
interaction, and the electric dipolar operator is so far taken for the whole ensemble of molecules, 
rather than as the dipole operator for the individual molecules. 



Figure 1. The ensemble of identical, similarly oriented, and mutually independent molecules. 

For many practical applications, however, the obtained general form is somewhat inconvenient 
when it comes to the numerical evaluation of the susceptibilities, since tables of wave functions, 
transition frequencies and their corresponding electric dipole moments etc. often are tabulated 
exclusively for the individual molecules themselves. Thus, we will now apply the obtained general 
theory to an ensemble of identical molecules, making a transition from the wave function and 
electric dipole operator of the whole ensemble to the wave function and electric dipole operator 
of the individual molecule, and find an explicit form of the electric susceptibilities in terms of the 
quantum mechanical matrix elements of the molecular dipole operator. 

We will now apply three assumptions of the ensemble of molecules, which introduce a significant 
simplification to the task of expressing the susceptibilities and macroscopic polarization density in 
terms of the individual molecular properties: 

1. The molecules of the ensemble are all identical. 

2. The molecules of the ensemble are mutually non- interacting. 

3. The molecules of the ensemble are identically oriented. 

From a wave functional approach, it is straightforward to make the transition from the de- 
scription of the ensemble down to the molecular level in a strict quantum mechanical perspective 
(as shown in Butcher and Cotters book). However, from a practical engineering point of view, 
the results are quite intuitively derived if we make the following observations, which immediately 
follow from the assumptions listed above: 

• The positions of the individual molecules does not affect the electric dipole moment of the 
whole ensemble in V, if we neglect the mutual interaction between the molecules. (This holds 
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only if the molecules are neutrally charged, since they otherwise could build up a total electric 
dipole moment of the ensemble.) 
• Since the mutual interaction between the molecules is neglected, we may just as well con- 
sider the individual molecules as constituting a set of "sub-ensembles" of the general ensemble 
picture. In this picture, all quantum mechanical expectation values, involved commutators, 
matrix elements, etc., should be considered for each subensemble instead, and the macroscopic 
polarization density will in this case become 

{the number of molecules 1 ( the expectation value of the 1 
within the volume V J \ molecular dipole operator J 

{the volume V} 





Figure 2. The general ensemble seen as an ensemble of identical mono- molecular sub-ensembles. 

Assuming there are M mutually non-interacting and similarly oriented molecules in the volume 
V, the macroscopic polarization density of the medium can be written as 



P(r,t) = ^(Q) = ^(-e^Tj+e^Zkh) 



V 



k 

nuclei 



electrons 

= {arrange terms as sum over the molecules of the ensemble} 

M 



^E(- e E f f )+e E^ m)f i ro) ) ; 



electrons 



clci 



molecular clcc. dipole operator, ef( m ) 
electric dipole moment of ensemble 



m—l 



(m)\ 



= {the molecules are mutually noninteracting} 
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{the molecules are identical and similarly oriented} 
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where N = M/V is the number of molecules per unit volume, and 

(ef) = Tr[g(t)er] 

is the expectation value of the mono-molecular electric dipole operator, with g(t) (that is to say, 
p(t) with a "kink" 1 to indicate the difference to the density operator of the general ensemble) is 
the molecular density operator. Notice that the form P(r,i) = N(er) is identical to the previous 
form for the general ensemble, though with the factor 1/V replaced by N — M/V, and with the 
electric dipole operator Q a of the ensemble replaced by the molecular dipole moment operator ef a . 

When making this transition, the condition that the molecules are mutually independent is 
simply a statement that we locally assume the superposition principle of the properties of the 
molecules (wave functions, electric dipole moments, etc.) to hold. 

Going in the limit of non-interacting molecules, we may picture the situation as in Fig. 2, 
with each molecule defining a sub-ensemble, which we are free to choose as our "small volume" of 
charged particles. As long as we do not make any claim to determine the exact individual positions 
of the charged particles together with their respective momentum (or any other pair of canonical 
variables which would violate the Heisenberg uncertainty relation), this is a perfectly valid picture, 
which provides the statistical expectation values of any observable property of the medium as M/V 
times the average statistical molecular observations. 

The previously described operators of a general ensemble of charged particles should in this case 
be replaced by their corresponding mono-molecular equivalents, as listed in the following table. 



General 


Molecular 




ensemble 


representation 






efft 


(Electric dipole operator) 


Po 


—> 9o 


(Density operator of thermal equilibrium) 


Pn(t) 


Qn(t) 


(nth order term of density operator) 


Ho 


— > Ho 


(Hamiltonian of thermal equilibrium) 


1/V 


TV = M/V 


(Number density of molecules) 


Po(a) 


Qo(a) 


(Molecular population density at state a}) 



Whenever an operator is expressed in the interaction picture, we should keep in mind that the 
corresponding time development operators Uo(t), which originally were expressed in terms of the 
thermal equilibrium Hamiltonian Ho of the ensemble, now should be expressed in terms of the 
mono-molecular thermal equilibrium Hamiltonian Ho- 

The matrix elements of the mono-molecular electric dipole operator, taken in the interaction 
picture, become 

(a\ef a (t)\b) = (a\U (-t)er a U (t)\b) 

= {the closure principle, B.&C. Eq. (3.11), [66'] ab = ^ O ak O kb } 

k 

= SS<o|^o(-*)|A;><A|ef a |/></|&o(t)|&> 

k I 

= {energy representation, B.&C. Eq. (4.54), [Uo(t)] ab = cxp(— iE a t/h)5 a b\ (3) 

= ^2 X! c *v{-i^a{-t)/h)5 ak {k\er a \l) exp(-iE;t//l)^& 
k i 

= (a|ef Q |6)exp[« (E„ - E b )t/h] 
' ^ ✓ 

= (a\ef a \b) exp(i£l ab t), 

1 kink n. 1. a sharp twist or bend in a wire, rope, hair, etc. [Collins Concise Dictionary, 
Harper-Collins (1995)]. 
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where fl a b = (E a — ¥.{,)/ h is the molecular transition frequency between the molecular states \a) and 
\b). We may thus, at least in some sense, interpret the interaction picture for the matrix elements 
of the molecular operators as an out-separation of the naturally occurring transition frequencies 
associated with a change from state \a) to \b) of the molecule. 
First order electric susceptibility 

From the general form of the first order (linear) electric susceptibility, as derived in the previous 
lecture, one obtains 

N f° 

X^2(-w;w) = -"r^ J Tr{£ [ef M ,er a (T)]}exp(-iwT)dT 



Ne 2 f° ^ 

J — oo „ 



eoih 

= {the closure principle} 
Ne 



e^iti 



I ^^(a|£o|fc)(fc|77^"( T ) ^ r Q (r)f M |a) exp(-i^T) dr 

J — oo „ 



k 

= {use that (a\g \k) = go(a)S ak } 



Ne 2 f° ^ 

= -— ^ J 2J ^o(a)[<o|^r Q (r)|a) - (a|f Q (r)f AI |a)] exp(-iwr) dr 



= {the closure principle again} 
Ne 2 ' 



E^°( a )E [( a M & )( & l^"( T )l a ) _ ( a l^«( r )l & )( & lvl a )] exp(-iwr) dr 

; a 6 

= {Use results of Eq. (3) and shorthand notation (a\f a \b) = r" b } 



Ne 2 
e^ih ■ 



E So{a) Y / [ r ab r ba exp(in ba r) - r a ah r^ a exp(-z0 6o r)] exp(-iwr) dr 

a b J - OG 

,0 

= {evaluate integral / exp(zCr) dr — > l/(iC)} 

J —CO 



Ne 2 ^ _ ^ ( ^ . r£ 6 r£, 



a 6 

This result, which was derived entirely under the assumption that all resonances of the medium 
are located far away from the angular frequency of the light, possesses a new type of symmetry, 
which can be seen if we make the interchange 

(-<*>, ^ 

which by using the above result for the linear susceptibility gives 

A/"„2 , a fJ- A 1 „a . 

= ^" E «>(«) E (rf^Z 2 - + 

CI o 



= {change order of appearance of the terms} 

a 



= X^(-w;w), 

which is a signature of the overall permutation symmetry that applies to nonresonant interactions. 
With this result, the reason for the peculiar notation with "(— u;;w)" should be all clear, namely 
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that it serves as an explicit way of notation of the overall permutation symmetry, whenever it 
applies. 

It should be noticed that while intrinsic permutation symmetry is a general property that 
solely applies to nonlinear optical interactions, the overall permutation symmetry applies to linear 
interactions as well. 

Second order electric susceptibility 

In similar to the linear electric susceptibility, starting from the general form of the second order 
(quadratic) electric susceptibility, one obtains 



xi 2 i(-^;wi,w 2 ) 

Ne 3 1 



s (ih) 2 2! 



Ne 3 1 



/ Tr{po[[^,r ct (Ti)],r / 3(r 2 )]}exp[-z(wiTi + uj 2 t 2 )] dr 2 dn 

-oo J — oo 
3 -| i-O /-Ti 

S / / E^ a ^°[ [^'^a( T i)] , r / 3(r 2 )]|a}exp[-z(wiTi + uj 2 t 2 )] dr 2 dr x 

J— oo J — oo „ ^ v ^ 



-oc ^ — oc 

r M r Q (T 1 )-r Q (T 1 )f M 



{expand the commutator} 
Ne 3 1 



S / / ^ Ua\Q r^r a {T{)rj3(T2)\a) - {a\Q Q f a (Ti)f^fp(T 2 )\a) 

J — oo J — oc a 

(r 2 )r M r Q (Ti)|a} + Hpo^/K^^tTi)^! )) exp[-i(wiri + w 2 t 2 )] dr 2 dri 



= {apply the closure principle and Eq. (3)} 
Ne 3 1 r° 



eo(ih) 2 2! 



S / I" EEEE (<a|eo|fc)(fc|v|6}<6|f Q (r 1 )|6>(c|f /3 (r 2 )|a) 

J — oo J —OO „ i, i „ x 



-oc^-oc a fc 6 c 

-( a |^o|fc}(fc|f Q (T 1 )|6}(6|f Al | C }(c|f /3 (r 2 )|a) 
- (a|^o|fc>(fc|r : /3(T2)|6)(6|f ([t |c)(c|f a (Ti)|a) 

+ (a|eo|fc}(fc|r/3(T2)|&)(fo|r Q (Ti)|c}(c|f M |a)^ exp[-i(a>iTi + w 2 t 2 )] dr 2 dri 

= {use that (a|g |^) = Qo{ a )S a k and apply Eq. (3)} 
Ne 3 1 — f° ^ 



e fiM2 ^j S E g °( a )EE / / | r a6 r fc>faexp[i(fi 6c ri + fi ca r 2 )] 

- Tabrb/L exp[i(Q a6 Ti + Q ca r 2 )] - r£ 6 r£.r£, exp[i(ft ca Ti + ft a6 r 2 )] 

+ r ab r bc r ca exp[i(0 6c Ti + O o6 T 2 )] | expf-^WiTi + W 2 T 2 )] dT 2 

]\f p 3 1 fO H a (3 

! „ V , ^-oo U (^ca - W 2 ) > ^ ' 



£o(^) 2 2! 

v ' a be 



exp [z (n ab + n ca ) n] - ff^ r ™ exp [i (o ca + o ab ) n )] 



i(O ca -W 2 ) v ' «(^o6-W2) s - 



ab be' ca 



+ - 



exp[i (Q bc + Q ab ) n] \ exp[-i (wi + w 2 ) n] dn 
* ' J * ' 



i(ft ab - uj 2 ) 

JVe 1„\-^ , ,\-^\-^f <„uluJ„„ 'ab'bc'ca 



eo( ^)2 2 , S E^( a )EE{ (n ca ~Zmba-UJ a ) + (fica-W2)(nc6-W„) 

r /3 r M r a r „a fj, - 
_| ab be ca ab' fee ca I 



(Q ab - Lo 2 )(n cb - uj a ) (Q ab - w 2 )(-O co - u a ) . 

. . . continued on next page . . . 
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. . . continuing from previous page . . . 



^£,o(a)££{; 



e H 2 2! ^ ( fiac + W 2)(^afc + W<r) (fi ac + w 2 )(£l hc + w CT ) 



n /5 „/3 r a r /i 
_| ab be ca 



Now it is easily seen that if we, for example, interchange 

^ (w 2 ,/3) 

one obtains 

(2) f \ 

!S£,o(a)££{ 



e h 2 2! ^ ^ . 1 -(^oc-w CT )(O a 6-W2) (n oc -w CT )(f2 6c -u; 2 ) 



„/3 „a H a fj - 
_ ab be ca | ab be' ca I 

(Q ba - LU a )(Qbc - W 2 ) (Cl b a - La a )(£l ca - Ld 2 ) > 

= {use Q ab = -Clba, etc.} 

PJ P 3 1 _ . . . . . r r P r ct r H r a J3 fi 

I\e 1 >— ^ ^ J 1 ab' be' ca ' ab' be' ca 



!S£^o(a)££{ 



£ fr 2 2! ^ u *■ (^ca + W CT )(f2 6a + W 2 ) (O co + W CT ) (fi ch + W 2 ) 



identify 4th term 



(i p a fJ- a 

r ab r bc r ca , r ab r bc 1 



(Q ab + uj a )(n cb + uj 2 ) (^ab + u a )(n ac + W 2 ) J 

s v ' 

identify 1st term 

= {interchange dummy indices a — > c — > 6 — > a in 2nd term} 
= {interchange dummy indices a ^ b ^ c ^ a in 3rd term} 

\T e 3 i , r P r ct r )i r a J3 ji 

Jve 1 « ' x^^^^l 1 ab' be' ca ' ca' ab' be 



^£,o(a)££{ 



e h 2 2! ^ u (^ca +w (T )(O ba + w 2 ) (£l fcc + W CT )(fi ba +CJ 2 ) 

a be „ v / v v ✓ 

identify 4th term identify as 3rd term 

r bc r ca r ab , r ab r bc r ca \ 



(tt bc + UJ a )(Q, ac + U) 2 ) {flab + W,j){Q. ac + U) 2 ) 



identify as 2nd term identify 1st term 

(2) < \ 

that is to say, the second order susceptibility is left invariant under any of the (2 + 1)! =6 possible 
pairwise permutations of (— {uj\,a), and (uj 2 ,(3); this is the overall permutation symmetry 
for the second order susceptibility, and applies whenever the interaction is moved far away from 
any resonance. 

We recapitulate that when deriving the form of the nonlinear susceptibilities that lead to 
intrinsic permutation symmetry, either in terms of polarization response functions in time domain 
or in terms of a mechanical spring model, nothing actually had to be stated regarding the nature 
of interaction. This is rather different from what we just obtained for the overall permutation 
symmetry, as being a signature of a nonresonant interaction between the light and matter, and 
this symmetry cannot (in contrary to the intrinsic permutation symmetry) be expressed unless the 
origin of interaction is considered. 2 



2 It should though be noticed that the derivation of the susceptibilities still may be performed 
within a mechanical spring model, as long as the resonance frequencies of the oscillator are removed 
far from the angular frequencies of the present light. 
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